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Bond-operator mean field equations for the square-lattice, S — 1/2 bilayer Heisenberg model 
are developed and solved numerically. In the vicinity of both the zero-field critical point and 
the field- induced transitions, comparisons are made with T — and finite-temperature strong 
coupling expansions. The mean-field theory suggests that the quantum critical region for the 
field-induced transitions is restricted to significantly lower temperatures than one might have 
concluded based on strong-coupling expansions or other numerical studies. 
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§1. Introduction 

The appearances of unconventional superconductivity and non-Fermi-liquid behavior near the 
quantum critical point of antiferromagnetic instability in strongly correlated electronic systems, 
such as copper oxides or heavy fermion systems, has attracted great interest in recent years. These 
observations have led to the investigation of the interplay between magnetic long-range order and 
novel quantum disordered states. Particular attention has been paid to low-dimensional quantum 
spin systems in which the classical ground states are destabilized by strong quantum fluctua- 
tionsJ3>>i>l)>i)|) which are experimentally realizable in several compounds. 

The bilayer Heisenberg model 

H = J ( E Sm-Su + £ s 2 , r s 2ij j + S M" S 2,* , (i-i) 

(taking the interlayer exchange as the unit of energy, and with running over nearest neighbor 
pairs on the square lattice) has attracted considerable interest&'B'&'l!''0'll3''l^l) as a testing ground 
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for notions of quantum criticality.E^'liil) O 



>n increasing J from zero, at T = the model passes from 



a spin-disordered ground state through a critical point to a Neel-ordered state — the same phase 
diagram as the nonlinear a model, with 1/J playing the role of the coupling parameter g. 

For S = 1/2, this model is directly relevant to the magnetic properties of several cuprate high-T c 
parent, antiferromagnetic compounds (and possibly the high-T c materials themselves©), and also 
to the spin-gap system BaCuSi206 

© 

If an external field 



is applied to the bilayer Heisenberg model, new T = phases and phase transitions are introduced. 
In particular, if at h = the system is in the disordered, gapped phase then there is a critical field 
h c , trivially related to the zero- field triplet gap A, since there is a simple Zeeman splitting of the 
triplet elementary excitations. For h > h c the gap is closed and the system exhibits algebraic order 
in the spin components transverse to the applied field. 

A great deal is known about this T = field-induced transition© In particular, d = 2 is the 
upper critical dimension, so critical exponents take on their classical values (modulo log corrections) 
in the bilayer Heisenberg model. In all dimensions above the lower critical dimension d = 1, the 
dynamic exponent z = 2 and the correlation length exponent v = 1/2. A variety of detailed 
model studies for one-dimensional spin- gap systems have been carried outJliP Work on two- 

dimensional systems exhibiting this transition is sparser,© and what is particularly not well known 
is how high in temperature the associated quantum critical region extends. This is a matter of 
some significance, if one is to ascribe the temperature dependence of properties of an experimental 
system to quantum critical behavior. 

In this paper, we study the field-induced transition in the bilayer Heisenberg model by means of 
bond-operator mean-field theory!) and strong-coupling expansions .0)'!!) Our principal conclusion 
is that the quantum critical regime is rather narrow, and that "quantum-critical-like" temperature 
dependences of quantities such as the magnetization are not sufficient to demonstrate that the 
system is in the true quantum critical regime. 

Bond operator mean-field theory has been applied to the spin-ladder (ID) and weakly coupled 
spin-ladder (3D) models® and shown to yield quantitatively correct results. Mean-field theory is 
a particularly appropriate tool for the problem at hand because the model is at its upper critical 
dimension, and so the asymptotic critical behavior will be correctly described up to logarithms. 
The strong-coupling expansions provide a good check of the accuracy of the mean-field theory, in 
the regimes where the former are reliable. 

In §2, we describe the bond-operator mean-field theory, where the bilayer Heisenberg model is 
described in terms of three types of interacting bosons with a local constraint. In §3, ground state 
phase diagram in the h-J plane and the critical properties of phase transitions are discussed. In 
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§4, thermodynamic properties of quantum disordered and quantum critical phases in a magnetic 
field are presented. The last section (§5) offers a summary and discussions. 

§2. Bond-Operator Mean-Field Theory 

Let us consider the quantum disordered phase of (1.1) in the limit J<1 which is the analytic 
continuation of the interlayer dimer Hamiltonian (J = 0). Then we describe the phase by intro- 
ducing bosonic degrees of freedom corresponding to local dimer states (singlet and triplet states) 
following the bond-operator representation of Sachdev and Bhatt:& 

| a >= s t| >=-L(|u>-UT>) (2-1) 
|^>=4|0>=--^(|TT>-|II>) (2.2) 

|ty>=4|0>=-^(|tt>+|U>) (2-3) 

|t,>=4|0>=-L(|U>+| IT>). (2.4) 

Here the four types of bosons satisfy bosonic commutation relations, and the left and right arrows 
in the kets represent the spin states on the first and second layers respectively. However there is 
a significant difference from usual bosonic systems because in the physical subspace only a single 
boson is allowed on each dimers. Thus the Hilbert space is restricted by a local constraint on the 
number operators of these bosons in each dimer 

4* + 4 A< = 1 (2-5) 

where the summation convention is used for Greek indices. 
In terms of these bosons, the spin operators in each layer can be expressed as 

1 

2 

Then the bond-operator expression of Hamiltonian (1.1) is obtained in the form 

3 1 

H = £(- 4*}* + 4*LA<) " E kw* + £A< - V 

i i 

~i E(! - M(4Aj<^ - ^A^A-i) • (2-8) 

Here we have introduced local chemical potentials m to account for the local constraint (2.5). 
Mean-field theory of the disordered phase is constructed by the following procedures. (1) An 



&2,i — s \ta,i ^a,i S i ^ e al3-y^/3^-y,i) ■ (2-7) 
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average chemical potential \x = /ij is introduced to approximate the local constraints by a global 
constraint. (2) Bose-Einstein condensation of singlet states is assumed: (s\) = (si) = s, where 
(• • •) denotes the thermodynamic average. (3) Terms containing four triplet operators are dropped. 
Here we would like to comment on the last numerically simplified assumption. The quadratic terms 
in ^-operators will in tern lead to a nonzero expectation value of {t a t a ). By taking quadratic 
decouplings of fourth-order terms in t a , we found that the inclusion of these terms changes the 
numerical results only slightly, even near the critical coupling (for example, this lowers J c discussed 
in the next section about 5%). 
The resulting mean-field Hamiltonian is easily diagonalized using the Bogoliubov transformation: 



ika = cosh 9 k t ka + sinh 8 k t_ ka (2.9) 

cosh 2 fe = -(Afc/wfc + 1) (2.10) 

HMFTi^s) =N(-^ S 2 - f IS 2 + ^-^)+J2 U; k(da^a + l) (2.11) 

k 

co k = v /A|-4A fc 2 (2.12) 

A fc = - - n + 2 Js 2 7fc , A fc = Js 2 lk (2.13) 

Ik = 7: (cos k x + cos ky) . (2-14) 



The only effect of including an external field is to Zeeman-split the excitations, so that instead 
of all three lying at uj k they have energies u k and u k ± h. 

The above mean-field description of the disordered phase will be valid if the temperature is much 
less than the minimum energy gap. However, at intermediate and higher temperatures the density 
of triplet bosons becomes sufficiently large that their interactions cannot be neglected. Troyer, 
Tsunetsugu and Wurtz@) have presented a simple prescription to modify the mean-field theory 
so that it correctly treats both the low and high-temperature limits. They noticed that the main 
problem of the bosonic description is its overcounting of entropy due to the global constraint on 
the number of elementary excitations. They reweighted the M-boson part in the partition function 
so that each multiplet contributes the correct entropy. In their formulation, the free energy of the 
triplet bosons in a uniform magnetic field is expressed (per dimer) as 

/ = ~ ln{l + [1 + 2 cosh(/%)M/3)} (2.15) 

where 

^) = ^E e ^ k • ( 2 - 16 ) 

k 

Using these formulae, one obtains self-consistent equations for the mean-field parameter s and 
chemical potential /i by minimizing the total free energy (F(s,fi) = /(s,/x) + constant(s, /i)) with 
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respect to those parameters: 



2 



3 



k 



VI + 2r 7fc 



3 + 



2(l + 2cosh(/3/t))e-^ k 
1 + [1 + 2cosh(/3/i)]z(/?) 







(2.17) 




2 2N ^ VI + 2r 7 fe 



3 + 



2(l + 2cosh(/3/i))e- /3u,k 
1 + [1 + 2cosh(/3/i)]z(/3) 







(2.18) 



where 



r = 2J S 2 /(4 " & ■ 



(2.19) 



It is straightforward to verify that this approach leads to the correct free energy and magnetization 
for isolated dimers (J = 0) for all temperatures and external fields. The numerical analysis of the 
above equations when J ^ is discussed below. 

§3. Ground State Phase Diagram 

In the classical (large-S") limit, the Heisenberg bilayer with h > has three types of ground 
states. One is the fully polarized ferromagnet, realized in the high-field region. The others are 
canted phases in which the z component of the magnetization is parallel to the external field to 
gain Zeeman energy, while x and y components of the spins are ordered to gain exchange energy. 
For quantum spins a disordered (dimer) phase appears around J = for sufficiently small h. As 
mentioned in Introduction, the phase diagram of (1.1) in zero magnetic field has been extensively 
studied. The order-disorder transition lies in the universality class of classical d = 3 Heisenberg 
model, and is associated with a dynamical exponent z = 1. However, the field-induced transitions© 
have dynamical exponent z = 2 and so should lie in the universality class of a d = 4 classical model. 
The reduction of symmetry from 0(3) to U(l) due to the magnetic field leads one to anticipate 
that the universality of the transition is d = 4 XY type. 

To estimate the phase boundary between dimer and canted phases, we performed T = strong- 
coupling expansions of the longitudinal susceptibilities transverse to the applied field up to the 8th 
order at the ordering vectors from interlayer dimer Hamiltonian (i.e. series expansions in powers 
of the intralayer couplings J about the interlayer dimer singlet state) using connected cluster 
method© By applying the differential approximant method0 assuming power-law divergence of 
the susceptibilities at critical points for fixed values of the applied field h, x ~ (Jc(h) — </)~ 7 , we 
obtained estimates for the critical lines J c (h) (and found associated critical exponents 7 close to 
the expected value of 1). We could also estimate phase boundaries using the triplet excitation gaps 
at zero field which had been studied by strong-coupling expansions!^ previously. These different 
estimates of the critical lines are consistent with each other, as shown in Fig. 1. The boundaries 
between the fully polarized phase and canted phases shown in that figure can be obtained from the 
lowest one magnon excitation gap in linear spin-wave theory. 
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Next we compare the results of the bond-operator mean-field theory at T = and h = with 
the results from the series expansions. At T = 0, the self-consistent equations (2.17), (2.18) for 
the two parameters s and \i can be reduced to a single equation for the parameter r introduced in 
(2.19), namely, 

r = J(5 - 3K (r)) (3.1) 

with 

K(r) = —Y ; 1 . (3.2) 

Once the parameter r is determined from the above set of equations, the mean-field parameter s 
and chemical potential fi are obtained from 

s *= h --l(E(r)+K{r)) (3.3) 

and 

H = -\ + ^- r {E{r)-K(r)) (3.4) 

where 

E ( r ) = jf E V 7 ! + 2r 7 fc • (3.5) 

As J grows from 0, the triplet gap at the wave vector k = (tt, tt) decreases, and the gap vanishes 
at J = J c , signaling the instability of dimer phase. The critical value J c can be obtained by setting 
r = 1/2 and J = J c in (3.1), giving the value J c = 0.437 which is somewhat larger than the 
series expansion estimate 0.393. In addition we obtain a critical spin-wave velocity 0.715 which is 
somewhat smaller than the series expansion estimate 0.744. The triplet gap for various values of 
J in the dimer phase, and excitation spectra for J = 0.2 and 0.35, were obtained from numerical 
solutions of the mean-field equations, with the results shown in Figs. 2 and 3. For negative J, the 
mean-field theory gives the same critical parameters, changing only the sign of J c . For comparison, 
the series expansion estimate for the negative critical coupling is© —0.433, so mean field theory is 
remarkably accurate in that case. 

In the following section, we restrict our discussion of several thermodynamic properties only for 
positive J but the qualitative results of thermodynamic behavior are the same for negative J. 

§4. Finite Temperature Properties 

Our principal purpose in this section is to discuss the finite temperature quantum critical prop- 
erties of the bilayer system, where a magnetic field is the tuning parameter and the system has a 
quantum disordered ground state at zero field. We analyzed these properties on the basis of the 
bond-operator mean-field theory as well as finite-temperature strong-coupling expansions up to 5th 
order (i.e. series expansions of the free energy in powers of the intralayer coupling J) using con- 
nected cluster method.© Calculations of the latter variety for some parameter cases have already 
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been carried out to 8th order by Elstner and Singh for bilayers both withEiP and withoutEJ an 
external field. It turns out that the results from directly summing 5th order series are not very 
different from those of 8th order series, at least in the regimes where we expect these calculations 
to be most reliable. 

Let us first discuss the thermodynamic properties of the system when h = 0. In Fig. 4, we 
show the energy gap A of triplet excitations as a function of temperature as obtained from the 
mean-field theory. For values of J less than J c , A is finite at T = and shows a sharp rise with 
increasing temperature and eventually tends to saturate, reflecting a rapid extinction of long-range 
correlations. However, at J = J c , A rises from zero linearly with temperature, consistent with the 
prediction of the non-linear a model, A ~ 1.O4T.0) We have calculated the specific heat Cy and 
uniform susceptibility x f° r J = 0-2 as functions of temperature, using both mean-field theory and 
series expansions, with the results shown in Fig. 5. Since the system has a quantum disordered 
ground state, the thermodynamic quantities show thermally activated behavior. For the specific 
heat, both methods yield nearly identical results. For the susceptibility, mean-field theory deviates 
by about 10% at intermediate temperatures from the results of series expansions. 

Let us now discuss the temperature dependence of properties of the system at a field-induced 
transition. In Figs. 6 and 7, we show the magnetization and specific heat as functions of T at 
the critical field h = h c = A for J = 0.2. In general, we expect power-law dependences of 
thermodynamic quantities at low temperatures if the parameters of a system are set to a quantum 
critical point. Mean-field theory gives a linear T dependence for both quantities in a very narrow 
low temperature region, T < 0.1 <C A. The agreement between two approaches is fairly good 
except for some artificial features seen in the summed series expansions at low temperatures. 

The series expansions at finite temperature should reduce to the strong-coupling expansion from 
the dimer singlet state in the limit T = 0. Thus the thermodynamic quantities obtained by direct 
summation of finite series necessarily decrease exponentially at sufficiently low temperature. How- 
ever we might have the possibility of seeing quantum critical behavior at intermediate temperatures 
if the tuning parameter is set equal or close to the quantum critical point of the infinite system. 
Thereby we must be careful to judge whether or not the system is really exhibiting quantum critical 
behavior in the relevant temperature regime. In Fig. 6, it is found that the magnetization starts 
to show apparently linear behavior towards T = at T ~ 0.3, but mean-field theory indicates that 
the system has not yet entered the quantum critical region, since at T = the magnetization has 
a different slope than in the intermediate-T "linear" regime. The specific heat behaves similarly, 
as seen in Fig. 7. (The series expansion result exhibits a pair of peaks and a dip but this should be 
regarded as an artifact due to low order of the expansion and the crudeness of the extrapolation 
technique.) 

The case we have presented above, for J = 0.2, is actually the one in which there is the best 
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chance of observing the true quantum critical regime for the field induced transition. As J — ► the 
quantum critical regime vanishes since there is no way that the critical properties will be exhibited 
for T > J. As J — > J c , h c — > and it is natural to expect that for T > h c (but not too large) the 
system will exhibit the quantum critical behavior associate with the zero-field transition. Thus the 
field-induced transition's quantum critical regime is squeezed at both ends. We have made rough 
estimates of the extent of that regime as a function of J, based on the behavior of M(T), with the 
results shown in Fig. 8. It therefore appears that observing the quantum critical regime for the 
field-induced transition will, in general, be very difficult to achieve either in numerical simulations 
or experiments. 

§5. Summary and Discussion 

We have investigated the ground state phase diagram and thermodynamic properties of bilayer 
Heisenberg model on the square-lattice in a uniform magnetic field. The ground state phase diagram 
of the model has been investigated by means of strong-coupling expansion around the dimer limit 
(considering both the gap and longitudinal susceptibility) and linear spin-wave theory. We found 
that the longitudinal susceptibility exponent for the field-induced transition appears to be very 
close to zv = 1, as expected. However, we cannot exclude the possibility of logarithmic corrections 
to scaling expected at upper critical dimensions, and further field-theoretical investigations© are 
needed to clarify this point. It might also be interesting to calculate the transverse susceptibility, 
for which one expects an exponent of {z + d)u = 2. 

Thermodynamic properties of the model were investigated also by means of bond-operator mean- 
field theory, and its results were compared to those of the finite temperature strong-coupling expan- 
sions. We found the theory provides a simple but reasonably good description of thermodynamic 
quantities at all temperatures. It appears that the quantum critical region for field-induced tran- 
sitions in this system is restricted to extremely low temperatures, even though the magnetization 
exhibits what might appear to be the expected quantum critical power law behavior at higher 
temperatures. 

Finally, we mention that a transition entirely analogous to the field induced transition in Heisen- 
berg bilayers can take place in bilayer quantum Hall systems© Recently, Troyer and Sachdev0) 
presented quantum Montecalro simulation on Heisenberg bilayers in a magnetic field to determine 
the universal Kosterlitz-Thouless transition temperature Tkt hi the vicinity of a zero-field quantum 
critical point. The universality of Tkt when the ground state is ordered leads us to an anticipation 
of the universality of a quantum critical crossover temperature when the ground state is disordered 
(h < h c = A) (which is not our estimated one but the one between field-induced quantum critical 
and quantum disordered regions). We leave this issue to a future study. 



8 



Acknowledgements 

This work has been supported by the U. S. National Science Foundation through grant DMR 
94-57928 (MPG). 

P. Fazekas and P. W. Anderson: Philos. Mag. 30 (1974) 423. 
M. P. Gelfand, R. R. P. Singh and D. A. Huse: Phys. Rev. B 40 (1989) 10801. 
S. Sachdev and R. Bhatt: Phys. Rev. B 41 (1990) 9323. 
N. Katoh and M. Imada: J. Phys. Soc. Jpn. 63 (1994) 4529. 
M. Troyer, H. Kontani and K. Ueda: Phys. Rev. Lett. 76, (1996) 3822. 
A. W. Sandvik and D. J. Scalapino: Phys. Rev. Lett. 72 (1994) 2777. 
K. Hida: J. Phys. Soc. Jpn. 61 (1991) 1013. 
M. P. Gelfand: Phys. Rev. B 53 (1996) 11309. 
Z. Weihong: Phys. Rev. B 55 (1997) 12267. 
T. Matsuda and K. Hida: J. Phys. Soc. Jpn. 59 (1990) 2223. 
K. K. Ng, F. C. Zhang and M. Ma: Phys. Rev. B 53 (1996) 12196. 
A. V. Chubukov and D. K. Morr: Phys. Rev. B 52 (1995) 3521. 
S. Chakravarty, B. I. Halperin and D. R. Nelson: Phys. Rev. B 39 (1989) 2344. 
A. V. Chubukov, S. Sachdev and J. Ye: Phys. Rev. B 49 (1994) 11919. 

A. J. Millis and H. Monien: Phys. Rev. Lett. 70 (1993) 2810; H. Monien and A. W. Sandvik: J. Low. Temp. 
Phys. 99 (1995) 343. 

Y. Sasago, K. Uchinokura, A. Zheludev and G. Shirane: Phys. Rev. B 55 (1997) 8357. 
S. Sachdev, T. Senthil and R. Shanker: Phys. Rev. B 50 (1994) 258. 
I. Affleck: Phys. Rev. B 41 (1990) 6697. 
T. Sakai: J. Phys. Soc. Jpn. 64 (1995) 251. 
K. Hida: J. Phys. Soc. Jpn. 67 (1998) 1416. 



N. Elstner and R. R. P. Singh fcond-mat/9803085| ) 



M. P. Gelfand, R. R. P. Singh and D. A. Huse: J. Stat. Phys. 59 (1990) 1093. 



N. Elstner and R. R. P. Singh ( :ond-mat/9710286) 



B. Normand and T. M. Rice: Phys. Rev. B 56 (1997) 8760. 

M. Troyer, H. Tsunetsugu and D. Wiirtz: Phys. Rev. B 50 (1994) 13515. 

M. E. Fisher and H. Au-Yang: J. Phys. A 12 (1979) 1677; D. L. Hunter and G. A. Baker, Jr.: Phys. Rev. B 
19 (1979) 3808. To carry out the calculations we have used the package DA written by Keith Briggs at the 
University of Adelaide. 

Y. Matsushita, M. P. Gelfand and C. Ishii: J. Phys. Soc. Jpn. 66 (1997) 3648. 
D. F. B. ten Haaf and J. M. J. van Leeuwen: Phys. Rev. B 46 (1992) 6313. 
S. Sachdev: Phys. Rev. B 55 (1997) 142. 



L. Zheng, S. Sachdev and S. D. Sarma: cond-mat/9709315 



M. Troyer and S. Sachdev: xmd-mat/9807393 



9 



FIGURE CAPTIONS 

Fig. 1. Phase diagram of bilayer Heisenberg model in the magnetic field. The canted phases 
are specified by the ordering vectors (7r,7r;7r) and (0,0; 7r), where the third component ir indicates 
antiferromagnetic orientation between layers. The open circles are series expansion estimates of 
the triplet excitation gaps-H* 

Fig. 2. Triplet gap as a function of J(> 0). The dots and solid line are series expansion esti- 
mates!^ and mean-field theory results, respectively. 

Fig. 3. Triplet excitation spectra for J = 0.2, 0.35. The dots and lines are series expansion 
estimates and mean- field theory results, respectively. 

Fig. 4. Triplet gaps for various J (< J c ) as a function of temperature T. 

Fig. 5. Uniform susceptibility and specific heat (per a dimer) as a function of temperature for 
J = 0.2 and h = 0. The dots and solid lines are the results of series expansion up to the 5-th order 
and mean-field theory, respectively. 

Fig. 6. Magnetization at the critical field for J = 0.2 as a function of temperature. The symbols 
are the same as in Fig. 5. The broken line is just a guide to the eye. 

Fig. 7. Specific heat at the critical field for J = 0.2 as a function of temperature. The symbols 
are the same as in Fig. 5. The broken line is just a guide to the eye. 

Fig. 8. Estimates of the extent of the quantum critical regime for the field- induced transition, as 
a function of J. The curved line is just a guide to the eye. 
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